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Why does the intermediate polar V405 Aurigae show a 
double-peaked spin pulse? 
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ABSTRACT 

V405 Aurigae is an intermediate polar showing a double-peaked pulsation in soft 
X-rays and a single-peaked pulsation in harder X-rays. From XMM-Newton observa- 
tions we find that the soft band is dominated by blackbody emission from the heated 
white-dwarf surface. Such emission is at a maximum when either magnetic pole points 
towards us, explaining the double-peaked pulsation. The symmetry of the pulses re- 
quires that the angle between the magnetic and spin axes be high. 

The single-peaked pulsation in harder X-rays is explained in the usual way, as 
a result of opacity in the accretion curtains. However, the high dipolc inclination 
means that the accretion curtains are nearly in the plane. Thus the outer regions of 
the curtains do not cross the line of sight to the accretion footprints, explaining the 
absence of the deep absorption dip characteristic of many intermediate polars. The 
sawtooth profile of this pulsation requires that the magnetic axis be offset from the 
white-dwarf centre. 

We remark also on the double-peaked optical emission in this star. We suggest 
that the difference between V405 Aur's spin pulse and those of other intermediate 
polars is the result of its high dipole inclination. 
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1 INTRODUCTION 

V405 Aurigae (RX J0558. 0+5353) was discovered in the 
Rosat All-Sky Survey and identified as an intermediate po- 
lar (a cataclysmic variable with a magnetic white-dwarf pri- 
mary) by Haberl et al. (1994). 

It is notable, firstly, for showing a soft blackbody com- 
ponent in the X-ray spectrum, one of a number of such 
objects discovered with Rosat. Secondly, its soft-X-ray and 
optical emission shows a double-peaked modulation at the 
white-dwarf spin period (e.g. Allan et al. 1996), whereas 
most of these stars show a single-peaked modulation (see, 
e.g., Patterson 1994 or Hellier 2001 for reviews of this class). 

The hard X-ray emission in intermediate polars (IPs) 
originates below a stand-off accretion shock near the mag- 
netic poles of the white dwarf. The soft blackbody emission 
is then understood as arising from heated white-dwarf sur- 
face around the accretion footprints. This is nearly always 
seen in the AM Her class of cataclysmic variable, but it is 
seen only in some IPs, for which the reason is unclear. 

The issue of why some IPs show a single-peaked pulsa- 
tion, whereas others show a double-peaked pulsation, is also 
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unclear. One idea (e.g. Hellier 1996; Allan et al. 1996; Nor- 
ton et al. 1999) notes that IPs with shorter spin periods will 
have smaller magnetospheres in which the accretion discs 
are disrupted nearer the white dwarf. This could result in 
shorter, fatter 'accretion curtains' of material which might 
have lower opacity in the vertical direction, thus preferen- 
tially beaming X-rays along magnetic field lines. The two 
magnetic poles would combine to produce a double-peaked 
pulsation. With longer spin periods, where disc disruption 
occurs further out, the opposite might hold, with tall, thin 
accretion curtains preferentially beaming X-rays out of the 
sides. The two poles would then act in phase, producing a 
single-peaked pulsation. 



The XMM-Newton X-ray satellite has a larger collect- 
ing area and better spectral resolution than Rosat, allowing 
us to return to V405 Aur with better X-ray data than previ- 
ously obtained. We report here on a 30-ks XMM-Newton ob- 
servation aimed at understanding the pulsation at the 545-s 
spin period of V405 Aur. 
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Figure 1. Power spectra of the X-ray data (MOS-1 + MOS-2) 
in the soft (upper panel) and hard (lower panel) energy bands. 



2 OBSERVATIONS AND LIGHTCURVES 

V405 Aur was observed by the XMM-Newton satellite 
(Jansen et al. 2001; Turner et al. 2001) on 2001 October 
5. A 30-ks observation was made with the EPIC-MOS cam- 
eras in timing uncompressed mode using the Thin Filter 
1; the PN camera was not in operation. The source was out- 
side the OM window, so these data were not used. 

We analysed the data using the xmm-SAS software 
v5.4.1. In timing mode the central CCD data is compressed 
onto a single axis, thus selection is based on columns rather 
than pixels. Only single or double pixel events with a zero 
quality flag were selected from a 44-column region centred 
on the source. We also removed column 315 from the MOS-2 
data, which showed a spurious signal at low energies. Owing 
to the limited coverage of the timing mode window we used 
an adjacent chip (CCD3) to estimate the background. 

The power spectrum of the 0.2-12 keV data from the 
MOS cameras is given in Fig.Q This shows prominent peaks 
at the spin frequency and its first harmonic (periods 545.5 s 
and 272.7 s respectively). We find that the first harmonic is 
the stronger at energies below 0.7 keV (hereafter the 'soft' 
band), while the fundamental is the stronger at energies 
above 0.7 keV (hereafter the 'hard' band), as previously re- 
ported by Allan et al. (1996). 

The folded lightcurves from the two bands (Fig. |5J 
confirm these results, showing that the soft emission has 
a double-peaked profile whereas the hard emission has a 
single-peaked, sawtooth modulation. 
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Figure 2. The X-ray data (MOS-1 + MOS-2) from the soft and 
hard bands (top panels), folded on the 545-s spin cycle. The lower 
panels contain the 0.2-0.7 / 0.7-12 keV and 0.7-4 / 5-12 keV soft- 
ness ratios. Typical error bars are shown. Phase zero is taken as 
the peak of the larger of the soft maxima, and corresponds to 
HJD 2452187.55904. 



3 SPECTROSCOPY 

Little timing mode data has yet been published, and we 
used pre-release versions of new quantum efficiency calibra- 
tion files and canned response matrices (Sembay, private 
communication). However, there were clear calibration dis- 
crepancies between the two EPIC-MOS instruments below 
0.4 keV (Fig. We thus added systematic errors of 8% to 
the data in this energy range. Other than this, we do not 
make allowance for calibration uncertainties in calculating 
X 2 values. The data from the two cameras were fitted simul- 
taneously in the fits reported below. 

The X-ray emission in an IP arises from plasma heated 
to X-ray temperatures at a stand-off accretion shock, which 
then cools as it approaches the white dwarf surface (e.g. 
Aizu 1973; Cropper et al. 1999). The CEMEKAL model re- 
produces a multi-temperature spectrum, with a power-law 
distribution of temperatures, based on the mekal plasma 
model. 

Fitting our phase-averaged spectrum with a CEMEKAL 
model plus simple and partial-covering absorption (as com- 
monly found in IPs) gave a poor fit (x^=20) owing to a large 
soft excess. Adding a blackbody emitter to model the soft 
emission improved the fit hugely to x^ = l-38. 

We then tried replacing the cemekal with two mekal 
components and found a significant improvement in fit qual- 
ity (x 2 =1276, x^=l-28) implying that a power-law is a poor 
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Component 


Parameter 
(Units) 


Value 


Error 


A bsn. 


nu (cm ) 


l.UO X 1U 


{ j-Ci no n 1 
^-t-u.uy, — u. ) 


Blackbody 


Kl ^KC V j 


Q (17 \y 1 C\ — 2 

o.yixiu 


{ J-Cl Qfi fl/lOA 

M~U. oo, — 




Normalisation 


1 34x 10 — 2 




A bsn 


1h (cm 2 ) 


9 fl£ v 1 f|21 


f J-f\ 9Q fl 90A 


Part. Absn. 


■Tin (cm" 2 ) 


6.08 xlO 22 


(+0.58, -0.52) 




CvrFract 


0.52 


(+0.03, -0.03) 


Gaussian 


Energy (keV) 


6.42 


(+0.02, -0.02) 




Normalisation 


1.96X10" 5 


(+0.44, -0.45) 




Eq. Width (eV) 


121 


(+27, -28) 


Mekal 


kT (keV) 


0.168 


(+0.008, -0.003) 




Abundance 


3.39X10- 2 


(+1.14, -0.62) 




Normalisation 


0.32 


(+0.16, -0.12) 


Mekal 


kT (keV) 


9.00 


(+0.92, -0.73) 




Abundance 


0.12 


(+0.35, -0.36) 




Normalisation 


1.31xl0~ 2 


(+0.08, -0.08) 



Table 1. Model components and parameters fitted to the phase- 
averaged spectrum (see Section 0. The errors are the 90% confi- 
dence errors according to formal statistics (to the same power of 
ten as the values); note however that these are likely underesti- 
mates, since they do not account for the calibration systematics. 



description of the temperature distribution in the plasma 
column. The best-fitting temperatures were 0.2 and 9 keV. 
The blackbody component was still necessary (xS=l-60 
without it), as were both simple and partial-covering absorp- 
tion (x?=l-71 without these). There were significant residu- 
als at the 6.4-keV iron fluorescence line, so we added a nar- 
row Gaussian component, reducing \ 2 to 1203 (xS=l-21). 

In the above models the simple and partial-covering ab- 
sorbers acted on all of the emission components. However, 
the partial-covering absorber is sufficiently dense as to com- 
pletely block any blackbody emission that it covers, thus its 
effect is redundant with a change in blackbody normalisa- 
tion. So as a final tweak we altered the model so that the 
blackbody was only absorbed by a simple absorber, which, 
for generality, was allowed to have a different column than 
the absorber acting on the mekals. This yielded a slight \ 2 
improvement (x 2 =H87, x£=l-19)- 

The resultant model is given in Tableland illustrated 
in Fig. 2] Note that the spectral parameters vary over spin 
phase, so those listed in Table^will be weighted averages. 



4 PHASE-RESOLVED SPECTROSCOPY 

We have found that below 0.7 keV the spectrum is domi- 
nated by a blackbody component, whereas above 0.7 keV 
the spectrum is dominated by optically thin plasma (Sec- 
tion |3J|. Further, the soft blackbody component is double- 
peaked over the spin cycle, whereas the higher-temperature 
plasma has a single-peaked, sawtooth pulse profile (Fig. 0. 
Thus, for phase- resolved spectroscopy, it makes sense to deal 
with the soft pulse and the hard pulse separately. 



4.1 The soft pulse 

To analyse the soft pulse we defined phase regions encom- 
passing the maxima (phases 0.41-0.61 and 0.89-1.09) and 




Energy (keV) 

Figure 3. The phase- averaged spectra from the MOS-1 and 
MOS-2 (circles) cameras. The systematic discrepancy below 0.4 
keV is clear. 



minima (phases 0.21-0.41 and 0.66-0.83), as illustrated in 
Fig.H 

We fitted all four phase regions simultaneously with 
the model developed above (Section [3J, which comprised 
a soft blackbody and two mekal emitters, with a simple 
absorber acting on the blackbody, and both simple and 
partial-covering absorbers acting on the mekal emitters. 
For analysing the soft pulse the fits were allowed to opti- 
mise for the entire spectrum, but the quoted \ 2 values were 
calculated only from 0.2-0.7 keV (the changes in \ 2 between 
models were similar if \ 2 was calculated over the entire en- 
ergy range). 

Allowing the model to optimise for each phase region 
separately yielded a best-fitting \ 2 of 263 (j£=1.27). Con- 
straining the simple absorber to be the same at all phases 
caused minimal change in fit quality but reduces the number 
of free parameters (x 2 =264, x 2 =1.25). 

Alternatively, allowing the absorption to vary freely, but 
forcing the blackbody emitter to remain constant between 
phase regions, yielded a x 2 of 308.51 (x 2 = l-44). This in- 
crease suggests that the soft modulation is a change in black- 
body normalisation rather than a change in absorption. 

Note that the 0.2-0.7 keV region is the worst affected 
by calibration uncertainties, so the quoted \ 2 values are un- 
reliable. However, if we fit the data from each camera sep- 
arately, without adding systematic uncertainties, we reach 
the same conclusion: that the soft pulse is mainly a change 
in visible emitting area, and not a change in absorption. 



4 Evans & Hellier 




l 10 

Energy (keV) 



Figure 4. The phase-averaged spectrum with the fitted model, and the contributions from each component. 



4.2 The Hard Pulse 

Since the 0.2-0.7 / 0.7-12 keV softness ratio is effectively 
a ratio of the blackbody emission to the plasma emission, 
we have also computed the 0.7-4 / 5-12 keV softness ratio, 
which is sensitive to spectral changes in the harder plasma 
emission. 

This ratio (Fig. shows little variation other than a 
hardening at phases ~ 0.0-0.2, near the maximum of the 
hard pulse. We thus extracted and compared spectra for the 
phase regions 0.0-0.2 and 0.4-0.9. 

Fitting the model to the two phase regions indepen- 
dently gave a total x 2 °f 1253 (x^=l-10; in this section 
X 2 values are for energies above 0.7 keV only; again, this 
does not affect the Ax 2 values). Constraining the absorp- 
tion (simple and partial covering) to have the same value in 
each phase region gave an identical x 2 of 1253, implying no 
absorption change. Alternatively, requiring the normalisa- 
tions of the mekals to be the same in each region worsened 
X 2 to 1530 { X l= 1.34). 



Thus we conclude that the hardening near flux maxi- 
mum is caused by a change in the ratio of the soft and hard 
mekal components. At all other phase regions there appears 
to be little or no change in the ratios of the two mekal com- 
ponents, nor in the absorption. This is despite the fact that 
the overall flux level changes considerably over these phase 
regions. Similar results, that the hard pulse does not show 
much energy dependence, were reported using BeppoSAX 
data by de Martino et al. (2004). 



5 DISCUSSION 

V405 Aurigae has two main spectral components in the X- 
ray band. Below 0.7 keV the spectrum is dominated by a 
blackbody component which shows a double-peaked modu- 
lation over the spin cycle. Above 0.7 keV it is dominated 
by thermal plasma emission which has a single-peaked, saw- 
tooth modulation over the spin cycle. 

IP X-ray modulations are often the result of absorp- 
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Figure 5. Schematic diagram showing the location of the emitting regions and magnetic poles at various phases. The dipole offset has 
been exaggerated. The arc-shaped regions are the blackbody emitting areas and the short lines depict the accreting field lines. 



tion as 'accretion curtains' of material sweep across the line 
of sight with the spin cycle (e.g. Hellier, Cropper & Mason 
1991; Kim & Beuermann 1995). Further, it has been pro- 
posed that whether the pulsation is single-peaked or double- 
peaked depends on whether the accretion column is short 
and fat, with greater horizontal opacity, or tall and thin, 
with greater vertical opacity (Hellier 1996; Allan et al. 1996; 
Norton et al. 1999). 

V405 Aur, however, does not fit the above model since, 
firstly, it shows both single-peaked and double-peaked be- 
havior in the same star. Secondly, our spectral modelling 
shows that the absorption varies little over spin phase (Sec- 
tion 0, which is the essence of the above model. 



5.1 Explaining the spin pulse 

We suggest that the modulation of the blackbody emission is 
simply the foreshortening of the accreting polecaps, which 
are viewed face on when a pole points towards us (maxi- 
mum) and foreshortened when the pole is near the white 
dwarf limb (minimum). Thus one of the poles is in the mid- 
dle of the visible face at each maximum (see Fig. 5). The 
difference between the two maxima is slight, so we require 
that both accretion polecaps be near the equator of the white 
dwarf, so that there is little difference in their visibility. This 
implies that the angle between the magnetic and spin axes 



is large, as is also suggested by the polarimetry reported by 
Shakhovskoj, Andronov & Kolesnikov (2001). 

Turning now to the hard X-ray pulse, the accretion ge- 
ometry for a highly inclined dipole is illustrated in Fig. |5] 
This shows that the difference in view between the accretion 
footprints at the upper and lower poles is relatively small. 
Nevertheless, we can adopt the standard accretion-curtain 
model (see, e.g., Hellier et al. 1991) and suppose that the 
lower pole viewed at phase is brighter than the upper pole 
viewed at phase 0.5 (Fig.|SJ. This is because the opacity in 
the column causes X-rays to emerge preferentially perpen- 
dicular to the curtain, which means that the lower pole is 
somewhat better presented. 

Note that with such a high dipole inclination, the cooler, 
outer regions of the accretion curtains never cross the line of 
sight. This explains why we don't see a prominent absorp- 
tion dip dominating the spin pulse, in contrast to many other 
IPs (e.g. AO Psc, Hellier et al. 1991). However, the X-rays 
will be emerging from a highly ionized post-shock region in 
which absorption and electron scattering will cause opac- 
ity, thus explaining the change in brightness with aspect of 
the accretion curtain, and the need for a partial-covering 
absorber in the spectral fit. 

Note that at phase the blackbody emission seen at 
the lower pole would be from heated surface on the side of 
the curtain further from the magnetic pole (that nearer the 
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pole being obscured by the curtain) whereas at phase 0.5 the 
blackbody emission at the upper pole is from the side nearer 
to the magnetic pole (see Fig. EJ. Thus an asymmetry be- 
tween these regions can explain the fact that one maximum 
is higher than the other. 



5.2 The sawtooth profile 

We further need to explain the fact that the hard- X-ray 
pulse shows a sawtooth profile, declining more rapidly than 
it rises. We suggest that this can be explained if the magnetic 
dipole is offset by ~ 0.1-0.2 white-dwarf radii. Such offsets 
are commonly found in magnetic white dwarfs (Wickramas- 
inghe & Ferrario 2000). 

From the deductions above, the lower pole is on the 
meridian facing us at phase 1 (thus presenting maximum 
blackbody area and producing a soft maximum). However, 
suppose that the dipole is offset to the left (Fig.|^J, so that 
the magnetic pole only points towards us shortly afterwards. 
This will slightly delay the hard X-ray peak. 

The lower pole then disappears over the limb at phase 
0.25. Because of the offset, the upper pole does not yet ap- 
pear, and this asymmetry explains the relatively rapid fall 
in flux at this phase. 

With the asymmetry pushing both poles to the far side 
of the white dwarf, we have lowest flux occurring before 0.5. 
Then the upper pole appears and is in the middle of the face 
near phase 0.5. Its less favorable aspect makes the system 
fainter than when we see the lower pole at phase 1.0. 

As we continue through the cycle towards phase 0.75, 
the offset of the dipole biases both poles towards the visible 
face of the white dwarf, which is why we see more flux ap- 
proaching phase 0.75 than we do following phase 0.25, thus 
explaining the sawtooth. However, since both accretion cur- 
tains are then presented with their azimuthal sweep being 
edge on to us, the flux level is lower than at phase 1.0. 



phase 0.5 phase 1.0 




Figure 6. Side-on schematic view of the model proposed in the 
text, shown here for a system inclination of 65° and a dipole 
inclination of 60°. 



5.4 The blackbody emitting area 

From our spectral fits we can estimate the blackbody emit- 
ting area. Note, however, that the calculated blackbody flux 
is strongly dependent upon the degree of absorption acting 
on it. In our best fitting model (Table^ this is surprisingly 
large (~ 10 21 cm -2 ), particularly if we associate the simple 
absorber with interstellar absorption. This in turn causes 
the flux measured to be high. However, forcing the column 
density to take on a lower value (e.g. 3xl0 20 cm -2 ) causes 
a worsening in x 2 °f <^ 100- Note, though, that the energy 
range over which the absorbed blackbody acts is that for 
which the calibration uncertainties are at their worst, mak- 
ing the values unreliable. In what follows we use the flux 
calculated from the best fit (rin = 10 21 cm -2 ). 

The blackbody emission component (Section [3J has an 
unabsorbed bolometric flux of 1.12xl0 -12 J m" 2 s" 1 with 
a temperature of 40 eV. This implies an emitting area of 
4.76x10 s (d/300 pc) 2 km 2 , which corresponds to ~8xl0~ 4 
of a white-dwarf surface [Haberl & Motch (1995) reported a 
corresponding value of 2x10 , based on a temperature of 
57 eV and a column of 5.7xl0 20 cm -2 .] This area is consis- 
tent with other estimates for the accretion area in IPs, for 
example the upper limit of 0.002 deduced by Hellier (1997) 
from eclipse timings of XY Ari. 



5.3 The optical spin pulse 

Although we do not present optical data in this paper, we 
can consider whether the model can explain the double- 
peaked nature of the optical spin pulse (e.g. Allan et al. 
1996; Skillman 1996; Still, Duck & Marsh 1998). 

In the standard accretion-curtain model, with a rela- 
tively low dipole inclination, the changing visibility of the 
curtains results in a single-peaked pulse with a maximum 
when the upper pole points away from us (Hellier et al. 
1991; Kim & Beuerman 1996). 

However, in V405 Aur we are postulating a large dipole 
inclination, such that the accretion curtains are nearly flat in 
the orbital plane (Fig.|SJ. Note that this is also supported by 
the optical polarimetry of Shakhovskoj et al. (2001) which 
implies equal visibility of two poles near the equator. 

With accretion curtains nearly in the plane, the view 
when the upper pole points towards us is similar to that 
when the lower pole points towards us (Fig. Thus we 
can explain a double-peaked optical pulse if the curtains are 
brighter when one of the poles points towards us and fainter 
when both are at quadrature (or vice versa). 



6 CONCLUSIONS 

XMM-Newton observations of V405 Aur confirms that the 
soft X-rays are dominated by a blackbody component, pre- 
sumably from heated white-dwarf surface surrounding the 
accretion regions. The heated region covers ~8xl0 -4 of 
the white dwarf, comparable to the upper limit found from 
eclipse timings of XY Ari. 

We propose that the modulation of this component re- 
sults from foreshortening of the blackbody emitting regions. 
We view the heated surface most favourably when either 
pole points towards us, so we see a double-peaked modula- 
tion in soft X-rays. The equality of the two maxima requires 
that the magnetic axis be highly inclined from the spin axis. 

The hard X-ray emission shows a single-peaked, saw- 
tooth modulation. This does not show the strong energy 
dependence of absorption, as usual in IPs. We suggest that 
this is because, with the high dipole inclination, the outer 
parts of the accretion curtains never cross the line of sight. 
However, electron scattering and opacity in the highly ion- 
ized post-shock column causes the intensity variation with 
spin phase. The sawtooth shape of the pulsation requires 
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that the magnetic axis be offset from the white-dwarf cen- 
tre. 

We also suggest that the high dipole inclination is re- 
sponsible for the double-peaked optical pulse. The high 
dipole inclination appears to be the main reason for the 
differences between V405 Aur's pulsation and typical IP be- 
haviour. 
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